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ABSTRACT 

In recent work we suggested that photons of energy >100 MeV detected from GRBs by the 
Fermi Satellite are produced via synchrotron emission in the external forward shock with a 
weak magnetic field - consistent with shock compressed upstream magnetic field of a few 
tens of micro-Gauss. Here we investigate whether electrons can be accelerated to energies 
such that they radiate synchrotron photons with energy up to about 10 GeV in this particular 
scenario. We do this using two methods: (i) we check if these electrons can be confined to the 
shock front; and (ii) we calculate radiative losses while they are being accelerated. We find 
that these electrons remain confined to the shock front, as long as the upstream magnetic field 
is > 10/iG, and don't suffer substantial radiative losses, the only condition required is that the 
external reverse shock emission be not too bright: peak flux less than 1 Jy in order to produce 
photons of 100 MeV, and less than ^100 mJy for producing 1-GeV photons. We also find 
that the acceleration time for electrons radiating at 100 MeV is a few seconds (in observer 
frame), and the acceleration time is somewhat longer for electrons radiating at a few GeV. 
This could explain the lack of >100 MeV photons for the first few seconds after the trigger 
time for long GRBs reported by the Fermi Satellite, and also the slight lag between photons 
of GeV and 100 MeV energies. We model the onset of the external forward shock light curve 
in this scenario and find it consistent with the sharp rise observed in the 100-MeV light curve 
of GRB080916C and similar bursts. 
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1 INTRODUCTION 

The Fermi Satellite has detected 18 GRBs (Gamma-ray Burst) at 
energies >100 MeV by LAT (Large Area Telescope). This emis- 
sion can be described as follows. The first 100 MeV photons ar- 
rive ~ 1 s (in the host galaxy frame) after the trigger time, for 
long GRBs; the trigger time is the time when low energy pho- 
tons (~ 1 MeV) are first detected by the GBM (Gamma-ray Burst 
Monitor aboard Fermi). The 100 MeV light curve rises fast until 
it peaks and then it decays as a single power-law for a long dura- 
tion of time (of order 10 3 s) - much longer than the duration of the 
lower energy photons detected by GBM - until it falls below the 
detector's sensitivity. Radiation above 100 MeV from GRBs has 
been suggested to be produced via the synchrotron mechanism in 
the external forward shock (Kumar & Barniol Duran 2009, 2010); 
the external forward shock scenario was first proposed by Rees & 
Meszaros (1992), Meszaros & Rees (1993), Paczynski & Rhoads 
(1993), and since then it has been used widely, see, e.g., Meszaros 
& Rees (1997), Sari, Piran, Narayan (1998), Dermer & Mitman 
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(1999), for a comprehensive review see, e.g., Piran (2004) and ref- 
erences therein. After our initial suggestion, many groups have also 
considered and provided evidence for this origin of the > 100 MeV 
radiation (Gao et al. 2009; Corsi, Guetta, Piro 2010; De Pasquale et 
al. 2010; Ghirlanda, Ghisellini, Nava 2010; Ghisellini, Ghirlanda, 
Nava 2010). The magnetic field required for this model is con- 
sistent with being produced via shock-compressed seed magnetic 
field in the CSM (circum-stellar medium) of strength of a few tens 
of micro-Gauss. The peak of the 100 MeV light curve can be at- 
tributed to the deceleration time which is the time it takes for the 
GRB-jet to transfer about half of its energy to the external medium. 

We investigate in this work whether electrons in the external 
forward shock can be accelerated to sufficiently high Lorentz fac- 
tors, even for a small CSM magnetic field of a few tens of /iG, so 
that the synchrotron radiation can extend to ~ 10 GeV as seen by 
FermifLKI for a number of GRBs. We study the electrons accel- 
eration in the context of diffusive shock acceleration (e.g., Krym- 
skii 1977, Axford, Leer & Skadron 1978, Bell 1978, Blandford & 
Ostriker 1978, Blandford & Eichler 1987), which was developed 
for non-relativistic shocks and has now been developed to consider 
relativistic shocks (semi-) analytically (e.g. Gallant & Achterberg 
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1999, Achterberg et al. 2001) and recently using 2D particle-in- 
cell simulations (e.g. Spitkovsky 2008a,b, Keshet et al. 2009). We 
assume that the electrons acceleration proceeds in the Bohm diffu- 
sion limit and that the magnetic field downstream is simply shock- 
compressed upstream magnetic field (other possibilities are con- 
sidered in, e.g., Milosavljevic & Nakar 2006, Sironi & Goodman 
2007, Goodman & MacFadyen 2008, Couch, Milosavljevic, Nakar 
2008). 

If the downstream magnetic field is simply the shock- 
compressed large-scale upstream field, then the field component 
perpendicular to the shock normal is amplified, while the paral- 
lel component is not. In this case, the downstream magnetic field 
will be mainly pointing to the direction perpendicular to the shock 
front normal, therefore particles trying to cross the shock front from 
downstream to upstream will find it difficult to catch up with the 
shock front, which moves with a speed of ~ c/3 with respect to the 
downstream medium (see, e.g., Achterberg et al. 2001, Lemoine, 
Pelletier & Revenu 2006, Pelletier, Lemoine & Marcowith 2009). 
One way that the particles might return to the upstream is if there 
is efficient cross-field diffusion of particles, which might occur if 
turbulent magnetic field is produced downstream (Jokipii 1987, 
Achterberg & Ball 1994, Achterberg et al. 2001). In principle, 
the turbulent magnetic field could dominate the shock-compressed 
field throughout the downstream region. However, it seems that al- 
though some turbulence is present just downstream of the shock 
front it does not persist across the entire downstream region (see re- 
cent simulations by Sironi & Spitkovsky 2010 that show that mag- 
netic field is amplified only right behind the shock front and re- 
turns to the shock-compressed value far downstream). In this case, 
much of the radiation is produced by particles swept downstream 
where the turbulence has died out and the magnetic field is consis- 
tent with the shock-compressed value. We also note that as long as 
the thickness of the turbulent magnetic field layer is smaller than 
the thickness of the shocked fluid divided by (Bt/Ba) 2 then the 
energy loss in the turbulent layer is small; B t is the turbulent mag- 
netic field strength and B d is the shock-compressed magnetic field. 
Therefore, in this work we neglect energy loss in the turbulent mag- 
netic field layer since it persists for a very short distance compared 
to the thickness of the shocked fluid (see, e.g., Keshet et al. 2009 
and references therein). 

This work is organized as follows. In Section 2 we address the 
question of high-energy electron confinement upstream and down- 
stream of the shock front, and also radiative losses suffered by elec- 
trons in between acceleration. Also, in Section 2, we discuss the lag 
of the >100 MeV light curves observed by Fermi LAT for several 
GRBs in light of our results on electron acceleration. In Section 3, 
we calculate the rise of the external forward shock light curve, tak- 
ing into consideration the non-zero time to accelerate electrons to 
high enough energies so they can radiate at > 100 MeV. We present 
our conclusions in Section 4. 



2 ELECTRON ACCELERATION FOR > 100 ME V 
EMISSION 

2.1 Electron confinement 

It is widely believed that electrons in non-relativistic shocks un- 
dergo diffusive shock acceleration, (e.g., Krymskii 1977, Axford et 
al. 1978, Bell 1978, Blandford & Ostriker 1978, Blandford & Eich- 
ler 1987). In the context of relativistic shocks, it has been shown 
that electrons gain energy each time they cross the shock front by a 



factor of ~ 2, except on the first crossing when they gain energy by 
a factor of the Lorentz Factor (LF) of the shock front (Achterberg 
et al. 2001). 

In order for electrons to turnaround while up/down stream and 
cross the shock front, their Larmor radius should be smaller than 
the size of the system, i.e. electrons should be confined to the sys- 
tem in order to be accelerated. In this subsection, we explore the 
confinement of electrons in the external forward shock model when 
the magnetic field in the unshocked medium, upstream of the shock 
front, is a few tens of fiG in strength, and the magnetic field in 
the shocked medium, downstream of the shock front, is simply the 
shock-compressed upstream field. 

The highest photon energy detected for Fermi GRBs is on 
the order of 10 GeV. We first calculate the random LF in the 
downstream co-moving frame, j e , of electrons radiating 10 GeV 
photons via synchrotron radiation, because these electrons have 
the largest Larmor radius and thus give us stricter confinement 
requirements. The synchrotron frequency in observer frame is 
Vs V n — eBd'f 2 T/2nm e c(l + z), where T is the bulk LF of the 
shocked fluid measured in the upstream rest frame (lab frame), 
Bd is the magnetic field downstream (measured in the local rest 
frame), z is the redshift, m e and e are the electron's mass and 
charge, respectively, and c is the speed of light (Rybicki & Light- 
man 1979). We convert the synchrotron frequency to 10 GeV, i.e., 
z-'io = hv S y n /l.G x 10~ 2 erg, where h is the Planck constant 
and 10 GeV corresponds to 1.6 x 10~ 2 erg. Using the convention 
Q x — Q/10 x and solving the last expression for j e yields 



7e = i.5xio 8 , 1 1 ( f(i + ^) 1/2 r 3 - 1 <- / 5 2 , 



(i) 



where B u is the magnetic field upstream, which is the magnetic 
field in the CSM. To obtain (1) we have assumed that the mag- 
netic field in the downstream region is Bd = 4TB U (Gallant & 
Achterberg 1999, Achterberg et al. 2001; note that the shock front 
LF measured in the lab frame is F s = V2F, Blandford & McKee 
1976), i.e. B d is the shock-compressed magnetic field in the up- 
stream (lab frame), which is what we have found for Fermi GRBs 
(Kumar & Barniol Duran 2009, 2010). 

The electrons' LF in the rest frame of the upstream plasma is 
7 e T, therefore, the Larmor radius in the upstream is given by 

2 

m e C 7el , OCw1n l9 x 1/2,. sl/2 D -3/2 

Rl,u = 5 = (2.6 x 10 cm) (1 + z) 1 B u _' 5 , (2) 

where we made use of (1) to eliminate f e . Comparing the Lar- 
mor radius with the size of the system upstream, R, which is 
given by the blast wave radius in the host galaxy rest frame — 
R = 2cr 2 t/(l + z) ~ 10 17 cm (where t ~ a few seconds and 
r ~ 10 3 is the blast wave Lorentz factor, e.g. Abdo et al. 2009a) 
— we find that Rl,u 3> R- This might suggest that electrons of 
7e ~ 10 8 are not confined to the system. However, an electron up- 
stream of the shock front travels only a distance ~ Rl,u/F before 
returning to the downstream, because by the time the angle between 
electrons' velocity vector and the normal to the shock front exceeds 
~ 1 /r, the shock front catches up with the electron and sweeps it 
back downstream (Achterberg et al. 2001). Therefore, for electron 
confinement upstream one should compare Rl,u/T with R: 



TR 



T 3 B. 



3/2 



1.1 



Ri 



Bl%(E 5i /n y/* : 



(3) 



where E is the energy in the blast wave, t is the time since the burst 
trigger in observer frame, and n is the number density of particles in 
the CSM; in deriving the second equality we made use of the time 
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dependence of Y and R in the external forward shock scenario for 
a homogeneous CSM (Sari, Piran, Narayan 1998). For B Ui -§ « 4 
found for the Fermi bursts, Rl,u/ (TR) < 0.2, and, thus, electrons 
radiating at 10 GeV cannot escape from the upstream side of the 
shock front; note that this conclusion holds for at least several hours 
in the observer frame. 

One should also check for electron confinement downstream. 
Here, the Larmor radius is smaller than it is upstream, because the 
magnetic field is larger by at least a factor of 4T due to shock com- 
pression. Therefore, the requirement for the confinement of elec- 
trons downstream is automatically satisfied whenever it is satisfied 
upstream. 

We conclude that there is no problem confining external for- 
ward shock electrons that radiate ~ 10 GeV synchrotron photons by 
the CSM magnetic field of strength > 10 fiG. 



2.2 Radiative losses during electron acceleration 

Electrons suffer radiative losses while being accelerated that could 
prevent them from reaching LFs of ~ 10 8 that are needed for ra- 
diating photons of 10 GeV via the synchrotron process. In this sec- 
tion, we ascertain whether or not the radiative losses suffered by 
electrons - due to synchrotron and inverse-Compton processes - 
are small compared with the energy gain in each round of crossing 
the shock front. We do this by comparing the total radiative cooling 
time-scale, t' cooh which is the time-scale for electrons to lose half 
of their energy, with the acceleration time-scale. 

For the case of ultra-relativistic shocks when the downstream 
magnetic field is simply the shock compressed upstream field, the 
upstream and downstream residency times for electrons are approx- 
imately equal, when particle diffusion is in the Bohm limit (Gal- 
lant & Achterberg 1999, Achterberg et al. 2001). Thus, the time 
it takes for electrons to make one complete cycle across the shock 
front is about twice the upstream residency time, and the upstream 
residency time is on the order of the gyro-time in the shock front 
co-moving frame (Baring 2004). In the lab frame, their upstream 
residency time is on the order of the time it takes them to travel 
a distance ~ Rl,u/Y. Since the Larmor radius (Rl,u) increases 
with increasing electron energy, the last shock crossing dominates 
the total upstream residency time. Thus, the time, in the co-moving 
frame of the blast wave, that electrons spend during the last cycle 
of crossing the shock front (upstream — > downstream — > upstream) 
before getting accelerated to Lorentz factor 7 e - given by (1) -is: 



~ = (1-7 x 10 3 s) + z^TfB^. 



(4) 



Taking into account the energy loss that these electrons expe- 
rience because of radiative cooling, the acceleration time-scale, in 
the blast wave co-moving frame, is given by 



*'acc(7e) ~ 4,(7e) + *' s (7e), 



(5) 



where ie 9 (7 e ) is the elapsed time since the beginning of the explo- 
sion when ig(7e) = t' CO oi(le)/2 (shock front crossing time should 
be equal to at least half of the the radiative cooling time in order to 
reach a particular y e ). At t' eq the electron barely reaches 7 e , there- 
fore, it needs an extra time on the order of ~ t' s to fully reach the 
desired j e . If t' s > t' cool /2, then the radiative cooling is too strong 
and prevents the electron from reaching the desired j e . In the sub- 
sections below we discuss synchrotron and inverse-Compton losses 
and calculate the radiative cooling time. 



2.2.1 Synchrotron losses 

The synchrotron cooling time-scale (in the blast wave co- 
moving frame) in the upstream of the shock front is t' syn u = 
6nm e c/aTB 2 -y e T 2 , where ar is the Thomson scattering cross- 
section. We find that the synchrotron cooling time for an external 
froward shock electron with LF given by (1) is 



(5.2 x io^r 1/2 (i+*r 1/2 r 3 - 1 iC / 5- 



(6) 



Since t' synyU > t' s by a factor of 30, then synchrotron cooling in 
the upstream is unimportant for electrons radiating at 10 GeV. 

Next, we calculate synchrotron losses in the downstream. 
Since t' 3yn oc B~ 2 , the synchrotron loss rate is larger downstream 
because of the larger magnetic field. For shock-compressed mag- 
netic field downstream, B is larger than upstream field by a fac- 
tor 4 (in the blast wave co-moving frame), and therefore t' synd = 
t' syn u /16. The effective synchrotron cooling time for electrons of 
LF given in (1) is t' syn « [l/2t' 3yn d + l/2t' 3y „ iU ] _1 , which gives 



t' syn = (6.1 x io 3 s) ^, 1/2 (i + z)- 1/2 r^B- 3 /l 



(7) 



We see from (4) that t' syn ~ 4t' a for electrons that produce 
synchrotron photons of 10 GeV energy, and therefore the maximum 
synchrotron photon energy — obtained by setting t' s = t' syn — is 
v m ax,syn ~ 40r3(l + z)^ 1 GeV (see, e.g., Guilbert, Fabian, Rees 
1983, de Jager et al. 1996, Cheng & Wei 1996). 



2.2.2 Inverse-Compton losses 

In this sub-section we calculate the inverse-Compton (IC) cooling 
time-scale for electrons. The inverse-Compton cooling time de- 
pends on the energy density of photons, and on the electron LF. 
Electrons in the external forward shock region are exposed to pho- 
tons from three different sources of radiation: (a) prompt ~MeV 
7-ray radiation which carries most of the energy release in GRBs; 
(b) synchrotron radiation produced in the external forward shock 
heated CSM and (c) radiation produced in the external reverse 
shock heated GRB-jet. We will consider all of these sources in our 
estimate for the IC cooling time. All calculations will be carried out 
in the rest frame of the shocked CSM. 
The IC cooling time is given by 

, _ 3m e c 2 

4 J dvaF'(v)^ £ 

where F'{y) is the energy flux in radiation per unit frequency in 
the co-moving frame of the shocked CSM, v is photon frequency 
in observer frame, and a is the cross-section for interaction be- 
tween electrons and photons; a « o~t (Thomson cross-section) 
when v < Fm e c 2 /[(I + z)hj e ] = Uk n , and for v 3> Vkn, 
a « erT^/^fcn) -1 . Thus, an approximate equationfor the IC cool- 
ing time is 



tic 



3m e c 

40T7e 



F'(< v kn ) + ^HF'(> v kn ) 



(9) 



where F'(< Vkn) is photon energy flux in the shock co-moving 
frame below the frequency v kn and F'{> v kn ) is the flux above 
Vkn- The frequency at the peak of the vF(y) spectrum is v v (in ob- 
server frame, i.e., co-moving synchrotron peak frequency boosted 
by a factor of V and redshift corrected) and Vkn, the Klein-Nishina 
frequency in the observer frame, for an electron of LF 7 e , is 
hvkn ~ (5eV) r 3 7 ( Tg(l + z) -1 . We note that for v kn > v p , only 
the first term in (9) should be kept. 
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The co-moving energy flux in radiation is related to the ob- 
served bolometric luminosity by: 



Lob 



Combining (9) and (10) we find 



t'ic 



3ttR T m e c 



1 + 



Vkr, 



(10) 



(11) 



where a is the spectral index, i.e. F'(y) oc v a for Vkn < v < v v \ 
for a > the term in the square bracket is of order unity. The 
above equation is valid only when v kn < v v . 

As mentioned before, there are three different sources of pho- 
tons that interact with electrons in the the external forward shock. 
We analyze these cases separately. 

Case (a): The prompt 7-ray emission in GRBs - the origin of which 
is still uncertain - often has a low energy spectral index a ~ 0, and 
the spectrum peaks at v v $ ~ 1. The luminosity of this component 
is the highest of the three cases considered; L (>s,52 ~ 10 for Fermi 
GRBs that have >10 2 MeV emission. The cooling time, obtained 
from (11), for this case is 



ic,. 



= (2.2xl0M- fe(1+Z) 



L b 



1 + 



(12) 



Case (b): The external forward shock synchrotron spectrum peaks 
at ~ lOOkeV (before the deceleration time), and the spectral index 
between vu n and v v is a ~ 1/3. The luminosity from the external 
forward shock is L oos ,52 ~ 0.1 at the deceleration radius (Rd), 
and at smaller radius it decreases as ~ R 3 . Therefore, we find from 
(11) that, for As? R d , 



t' IC , b w (2.2 x 10 4 s)i? 17 1 r 3 ^, 5 i ob 1 s 51 Rl 17 (l + z) 



(13) 



Case (c): If the GRB-jet is composed of protons and electrons, 
then the interaction of the jet with the CSM will heat up these par- 
ticles by the reverse shock propagating into the cold jet, and the 
synchrotron radiation produced would be very effective at cooling 
electrons in the forward shock region. This is because the peak of 
the reverse shock emission at the deceleration time is typically at a 
few eV (Sari & Piran 1999a), which is of order v kn for electrons 
of 7 e ~ 10 s . Since v v ~ v k n, then we can keep only the first 
term in (9), and use (10) for flux in the calculation of the cooling 
time. The observed luminosity (at the deceleration time) is given by 
L bs,d ~ 4ird%i'p l dFp l d, where dh is the luminosity distance, and 
F p ^d and v Pt d are the observed external reverse shock flux and peak 
energy at the deceleration time, respectively. Thus, the IC cooling 
time-scale for electrons in the external forward shock region, due 
to the radiation produced in the reverse shock heated GRB-jet, is 



t'ic,c w (400s)< 17 r| 3 dZ? 28 7e~>;;d^ 



-1 

,<2> 



(14) 



where Td is the LF of the GRB-jet at the deceleration time, F Pi d is 
in Jansky (Jy) and v p ,d in eV; the reverse peak flux can be ~ 1 Jy 
for very bright bursts such as GRB 990123 (Sari & Piran 1999b). 
The total IC cooling time is 



t', 



ic 



1 



t'r, 



t'r. 



t'r. 



(15) 



7C,« v IC,b IC,c_ 

and, finally, the total radiative cooling time, t' cooh is given by 
-1 



t'r. 



cool 



ff 



t'r 



IC 



(16) 



GRB 




v \MeV\ 


D ii,-5 


jjprom.pt 
obs, 53 


L ES 
obs, 51 


0809 16C 





0.5 


4 


1.5 


30 


090510 


0.4 


2.8 


2 


3.6 


1.6 


090902B 


0.4 


0.7 


2 


1.2 


1.7 



Table 1. The main quantities used in our analysis for three Fermi GRBs. a 
is the approximate spectral energy index, during the prompt emission phase, 
below the peak of the spectrum (/„ oc v" for u < v v ), v v is the observed 
peak of the spectrum; B„ -5 is the average upstream magnetic field, in 
units of 10 ^iGauss, obtained by modeling of the data for these bursts (Ku- 
mar & Barniol Duran 2010); L v J b ° r ^T and L^ b s g ri are the approximate 
observed isotropic equivalent luminosities of the prompt 7-rays and exter- 
nal forward shock emission at the deceleration time, respectively. Data are 
taken from Abdo et al. (20()9a,b,c). B u was obtained by setting three simple 
constraints while modeling the external forward shock emission: 1. Its flux 
at 100 MeV should agree with the observed value, 2. Its cooling frequency 
should be below 100 MeV at early times for consistency with the observed 
spectrum, and 3. Its flux at 100 keV should be smaller than the observed 
value during the observed steep decay, so that the external forward shock 
emission doesn't prevent the 100 keV to decay steeply (Kumar & Barniol 
Duran 2009, 2010). 



2.3 Application to Fermi GRBs 

In this section, we analyze 3 GRBs detected by Fermi: 
GRB080916C, 090510 and 090902B (Abdo et al. 2009a,b,c). The 
relevant data for each burst are tabulated in Table 1. We apply the 
above general results to these three GRBs, and determine the time 
it would take for electrons in the external forward shock for these 
bursts to be accelerated - via shock acceleration - to LFs capable 
of producing synchrotron photons of energies 100 MeV and 1 GeV 
(Table 2). 

The external reverse shock emission depends on the highly 
uncertain magnetic field strength in the GRB-jet, and it is there- 
fore difficult to estimate with any confidence. We calculate t' cool by 
neglecting the contribution of reverse shock emission to inverse- 
Compton cooling of electrons (t' ICc ), and this provides a lower 
bound to t' acc which is reported in Table 2 as a fraction of the de- 
celeration time, t' d — (1.7 x 10 3 s) Rd,i7^d3' f° r severa l Fermi 
bursts. We also provide in Table 2 an upper limit for the exter- 
nal reverse shock peak flux that is obtained by the condition that 
t'cool = t' s , at the deceleration time, when the contribution of the 
external reverse shock emission is included in the calculation of 

L cool- 

GRB080916C: The first MOO MeV photons arrived ~ 3s after the 
trigger time and then the 100 MeV light curve rose rapidly, as ~ t 6 , 
and peaked at ~ 5s (Abdo et al. 2009a). After the peak, which we 
identify as the deceleration time, t d , the flux decayed as a single 
power-law (this power-law is consistent with the expectation of the 
external forward shock model). So the first >100 MeV photons 
arrived att/td ~ 0.6, and photons of energies >GeV were detected 
at ~ 7 s (t/t d ~ 1.4). The highest energy photon, ~ 13 GeV, was 
detected ~ 16 s after the trigger time (t/td ~ 3). 

For electrons to produce 100 MeV photons their LF should 
be ~ 10 7 for this burst, and for 1 GeV photons the required 
7e > 3 x 10 7 ; we used B Uj _5 ~ 4 as suggested by the data for 
this burst (Kumar & Barniol Duran 2009) - see Table 1. The ac- 
celeration time for electrons to attain these LFs is calculated using 
(5); note that our theoretical estimates are roughly consistent with 
the observed time-scales for GRB080916C to within a factor ~ 2 
uncertainty of our estimates (Table 2). 
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Expected Observed 



GRB 


V 


tacc/td 


t/t d 


j?max 


osoyioc 


100 MeV 


0.3 


0.6 


0.30 




1 GeV 


0.6 


1.4 


0.02 


0905 10 


100 MeV 


0.3 


0.5 


9.90 




1 GeV 


0.6 


1 


0.90 


090902B 


100 MeV 


0.3 


0.3 


1.20 




1 GeV 


0.7 


1 


0.10 



Table 2. t aC c/td is the ratio of the time for electron acceleration to a spe- 
cific energy (corresponding to synchrotron frequency given in column 2) 
and the deceleration time; it is a measure of the delay, with respect of the 
trigger time, for photons of a given energy to arrive at the observer when 
the external reverse shock emission is smaller than given in the last column 
of the Table. The observed time delay of photons in column 2 with respect 
to the trigger time is t/td- F m £ x (in Jy) is the maximum possible observed 
external reverse shock peak flux, so that electrons can be accelerated to pro- 
duce photons of energy given in column 2alt d - 



GRB090510: For GRB0905 10 (Abdo et al. 2009b) there was a short 
delay in the detection of >100 MeV photons by ~ 0.1 s (we take 
the trigger time to be ~ 0.5 s after the GBM trigger, because of the 
presence of a precursor). The 100 MeV light curve peaked at ~ 0.2 
s (which we associate with the deceleration time), and so the arrival 
of the first >100 MeV photons was at t/td ~ 0.5. Higher energy 
photons arrived later: >1 GeV photons started arriving at td, and 
~ 10 GeV photons arrived slightly after t d - As shown in Table 
2 these results are roughly consistent with our estimates within a 
factor of 2. 

GRB090902B: The 100 MeV light curve for this burst peaked at 
~ 10s, which we identify as td, and the first >100 MeV photons 
were detected at ~ 3s after the trigger time (Abdo et al. 2009c), 
i.e. t/td ~ 0.3. Most of the GeV photons arrived at ~ td- The first 
10 GeV photon is detected at ~ 12 s. The highest energy photon 
detected was ~ 30 GeV at 80 s, i.e. at ~ 8tdEl The arrival time for 
the first > 100 MeV photons from this burst agrees with the electron 
acceleration time (Table 2). 

To summarize the main results of this section, it takes a few 
seconds for electrons in the external forward shock to be acceler- 
ated to a LF so that they can produce 100 MeV photons, and it takes 
a bit longer time for them to produce GeV photons. For this reason, 
GeV photons lag the 100 MeV radiation. If the external reverse 
shock flux is high (~ 1 Jy), then the first 100 MeV photons will 
be detected after the deceleration time, and 10 GeV photons will 
be detected much later (~ I0t d ), when the reverse shock flux has 
decreased substantially . If the external reverse shock flux is small 
(~ 10 mJy), then the first 100 MeV photons will arrive at about 
a third of the deceleration time, and GeV photons will be detected 
starting from close to the deceleration time. 



1 At this time, the LF has dropped by a factor of 8 3 / 8 ~ 2, and at z = 1.8, 
Vmax,syn ~ 10 GeV, a factor of ~ 4 smaller than the observed value. It 
can be shown that inhomogeneous magnetic fields lead to an increase of 
Vmax.syn by about an order of magnitude. 
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Figure 1. Top. The expected external forward shock light curve when 
the non-zero acceleration time of the emitting electrons is taken into ac- 
count. We plot the specific flux (normalized to the flux at td, f d ) versus 
time (normalized to the deceleration time, td)- We show two cases: (1) 
Ro < Rd < Rf (red solid line), and (2) Rd < Ro < R f (blue dashed 
line). Bottom. The light curve temporal slope, f3 = d ln(f / fd)/dln(t/td). 
The horizontal black dotted line shows the asymptotic value of the temporal 
decay index if we take electrons to accelerate instantaneously; for t < td 
the light curve would rise as i 3 . However, we find that the external forward 
shock light curve rises faster than t 3 due to the finite time it takes electrons 
to accelerate. 



3 STEEP RISE OF THE HIGH-ENERGY PHOTON 
LIGHT CURVE 

In this section, we calculate the onset of light curves of high-energy 
photons (>100 MeV). According to the standard external forward 
shock model, and assuming instantaneous acceleration of electrons, 
the observed flux rises as t s when the CSM-density is homoge- 
neous. We show here that the light curve rises much more steeply 
- similar to what is seen by FermifLAT data - when finite time for 
electron acceleration is taken into consideration, as was suggested 
by Kumar & Barniol Duran (2009). 

We calculate the rise of the light curve using a simple model. 
The external shock emission at some frequency v is zero until the 
blast wave reaches a radius Ro, which is set by the time-scale for 
electrons to be accelerated to a LF so that they start radiating at v; 
this time is calculated in Section 2. Electron distribution function 
in the neighborhood of the desired LF is assumed to grow with 
radius as oc R x , and the distribution attains its asymptotic power- 
law shape at some radius Rf ~ 2Rq. The rise of the light curve 
depends on Ro, Rf, x and the deceleration radius; the rise also 
depends weakly on the type of CSM and the energy spectral index. 

We show in Fig. 1 light curves for two different regimes Ro < 
Rd < Rf (case 1) and R d < Ro < Rf (case 2); for Ro < Rf < 
Rd the light curve is similar to case 1 except that between t f (the 
observer time corresponding to radius Rf) and td the light curve 
rises as ~ t 3 . Guided by the estimates provided in §2 we take Ro ~ 
R d /3 for 100 MeV photons, whereas for 10 GeV R ~ R d - We 
show the results for these choices of parameters in Fig. 1; note the 
very steep rise of light curves which appear similar to the fast rise 
of the observed 100 MeV light curve for GRB080916C as reported 
in Abdo et al. (2009a). 
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4 CONCLUSIONS 

In this paper we have investigated the acceleration of electrons 
via diffusion shock acceleration in the external forward shock of 
GRBs, and its implications for the high-energy photon detection by 
the Fermi Satellite. The external shock model, with a weak mag- 
netic field, has been proposed as the origin of the observed >100 
MeV emission detected by the Fermi Satellite from a number of 
GRBs (Kumar & Barniol Duran 2009, 2010). We find that high- 
energy electrons of Lorentz factor ~ 10 8 , required for producing 
-~10 GeV photons via the synchrotron process, can indeed be ac- 
celerated in an external shock that is moving through a CSM with a 
magnetic field of strength a few tens of fiG; they remain confined to 
the shock front as long as the upstream magnetic field is > 10/iG. 

We have also calculated the time it takes for electrons to be 
accelerated to a Lorentz factor ~ 10 7 so that they can radiate syn- 
chrotron photons at ~100 MeV. We find this acceleration time to 
be a few seconds in the observer frame; this calculation took into 
account radiation losses suffered during the acceleration process. 
This result offers a straightforward explanation as to why, for most 
Fermi GRBs, 100 MeV photons are not observed right at the trigger 
time but a little later. This also explains, why 100 MeV photons are 
observed before GeV radiation: it takes electrons radiating at GeV 
energies even longer time to accelerate. Taking this acceleration 
time into consideration while calculating high-energy light curves, 
we find that the light curve rises very rapidly - much faster than 
it does for the external forward shock model with instantaneous 
electron acceleration for which the flux rises as i 3 when the CSM 
has uniform density (the t 3 rise reflects the increasing number of 
swept-up electrons before the blast wave decelerates). 

The detection of the first 100 MeV photons at some fraction of 
the deceleration time, the longer delays in the detection of higher 
energy photonjf] and the fast rise of the 100 MeV light curve, fol- 
low the expectation of the external forward shock model when the 
finite time for electron acceleration is taken into account. Detec- 
tion of synchrotron photons of different energies provides an upper 
limit for the radiation flux produced in the reverse shock heated 
GRB-jet. For instance, the peak flux for the external reverse shock 
emission — if the peak of the spectrum is at a few eV — couldn't 
have been larger than about 300 mjy close to the deceleration time, 
for GRB080916C, otherwise it would prevent electrons from ac- 
celerating to a Lorentz factor of ~ 10 7 so that they can produce 
synchrotron photons of 100 MeV energy at early times (see Table 
2). Similarly, the reverse-shock flux should be < 20 mjy for GRB 
0809 16C in order that electrons in the forward shock are acceler- 
ated to a LF so that they produce 1 GeV photons. 

We speculate that the lack of >100MeV emission during 
the prompt phase of GRBs might be due to the presence of a 
bright optical source with observed flux larger than about 100 mjy, 
which would prevent electrons from reaching high Lorentz factors. 
This, coupled with the fact that GRBs with the largest LFs, which 
have small deceleration time, are the most likely bursts to be de- 
tected by Fermi (Kumar & Barniol Duran 2009) might explain the 
detection/non-detection of > lOOMeV radiation from GRBs. 

We note that the shock-compressed magnetic field scenario 
requires some cross-field diffusion of particles - presumably gen- 
erated by turbulence - to allow them to travel back to the upstream 

2 Note that this possible trend in the data goes in the opposite direction 
than in the prompt ~ 1 MeV emission, where higher energy photons arrive 
earlier than lower energy photons in long GRBs and there is no lag detected 
for short GRBs (Norris et al. 1986, Norris & Bonnell 2006). 



(e.g. Achterberg et al. 2001, Lemoine et al. 2006). This turbulent 
layer probably occupies a small fraction of the downstream region 
as suggested by recent simulations by Sironi & Spitkovsky (2010). 
Therefore, the picture that seems to emerge from numerical simula- 
tions and Fermi observations, is that there might be a small region 
of turbulence behind the shock front that aids in the acceleration 
of particles across the shock, but that the radiation is mainly pro- 
duced by particles that are swept downstream where the value of 
the downstream field is consistent with simple shock-compression 
of upstream field. 

There exists also the possibility that the CSM seed field is ac- 
tually a few /iG and some instability produced ahead of the shock 
amplifies it to the value of a few tens of /iG we infer by our mod- 
eling of Fermi GRBs (Kumar & Barniol Duran 2009, 2010). These 
instabilities have been studied by, e.g., Milosavljevic & Nakar 
2006, Sironi & Goodman 2007, Goodman & MacFadyen 2008, 
Couch et al. 2008. However, this possible amplification of a factor 
of ~10 is much smaller than the amplification customarily invoked 
to explain afterglow observations. 

We received a preprint from Piran & Nakar (2010) soon after 
this paper was completed. They have also considered the accelera- 
tion of electrons in the external shock. 
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